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Abstract. Combining the relative vicinity of the Local Group spiral galaxy 
CN1 ■ M 33 with the Spitzer images, we investigate the properties of infrared (IR) emis- 

q ' sion sites and assess the reliability of the IR emission as a star formation tracer. 

, We compared the photometric results for several samples of three known types 

^"*^ ■ of discrete sources (Hll regions, supernovae remnants and planetary nebulae) 

with theoretical diagnostic diagrams, and derived the spectral energy distribu- 
tion (from 3.6 to 24 fim) of each type of object. Moreover, we generated a 
catalogue of 24 /im sources and inferred their nature from the observed and 
^ c"| theoretical colours of the known type sources. We estimated the star formation 

O |. rate in M33 both globally and locally, from the IR emission and from the Ha 

emission line. 
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1. The IR emission of nebulae in M 33 
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CNJ ■ In the four IRAC and 24 /im MIPS bands, we performed photometry of cata- 

logued discrete sources of various type: Hn regions, SNRs and PNe. To better 
understand the IRAC/MIPS colour diagrams of discrete sources, we assembled 
a series of theoretical models (blackbodies, SB99, Cloudy, DUSTY. A PAH- 
dominated spectrum such as the diffuse ISM and the small dense Cloudy Hn 
region have virtually identical IRAC colours to some of the T n =300 K DUSTY 
5_i ■ models. The implication is that strong PAH features can be easily confused 

with a rising continuum in the MIR. The only way to resolve the ambiguity is 
through detailed spectral information or by using the 24 fJ-m filter in conjunction 
with the IRAC 8 one. Indeed, only 2 IRAC filters, together with MIPS-24, 
are sufficient to separate strong continuum emission in the MIR from strong 
PAH features. 

In the colour diagram of the known type sources, we observe a continuity 
among the various types: the Hn regions appear more dusty than the PNe, 
which behave more like stellar objects. Hll regions are better modelled by pure 
Cloudy Hn regions (including PAHs) and by DUSTY models with inner-shell 
temperatures around 600 K and low extinction. Statistically, the Hn regions are 
more concentrated near the regions defined by low values of the two shortest 
IRAC channels and/or equivalently high values of the two longer IRAC wave- 
lengths. The most intense PAH bands (6.2, 7.7, 8.6 /xm) are located in the 5.8 
and 8.0 /mi channels and this may explain the colours of Hn regions. 
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2. The 24 fim emission and star formation 

Spitzer images provide a unique opportunity to compile vast catalogues of MIR 
discrete sources down to stellar luminosities for galaxies in the Local Group: we 
extracted 515 sources at 24 /im (24Ss). In the colour diagrams, the region with 
the highest concentration of 24Ss matches well the one of Hn regions, although 
some of the 24Ss extend into the regions populated by SNRs and PNe. However, 
neither a pure blackbody, Cloudy, or the two DUSTY models can alone repro- 
duce the colour distribution of the 24Ss. The cumulative luminosity function of 
the 24Ss displays a double slope behaviour, markedly steeper at the high lumi- 
nosity tail. The change of slope, in the simplest scenarios, represents the change 
of regime between poor and rich clusters. The luminosity distribution appears 
to precisely match the whole range expected for ionised regions and their com- 
plexes: from the faintest B2 stars (0.1 mJy) to the brightest clusters containing 
more than 1000 ionising stars (1.8 Jy). We compare the star formation rates in- 
ferred from the MIR emission with that derived from Ha emission (uncorrected 
from extinction). At large scale, the two methods give a SFR of the inner disk 
(within 5 kpc) of about 0.2 M Q yr _1 . The relation between 24 /jm emission and 
Hll regions is also investigated by searching directly for Ha counterparts of the 
24Ss: the SFR for the single sources in IR and Ha are consistent but the scatter 
is rather large and not imputable to extinction alone. 

3. Conclusions 

The colours of the typical IR emissions of Hn regions, supernovae remnants and 
planetary nebulae are continuous among the different samples, with overlapping 
regions in the diagnostic diagrams. The comparison between the model results 
and the colours of Hn regions indicates a dusty envelope at relatively high tem- 
peratures ~600 K, and moderate extinction Ay< 10. The 24 /im sources IR 
colours follow the regions observationally defined by the three classes of known 
objects but the majority of them represent Hn regions. The derived total IR lu- 
minosity function is in fact very similar to the Hn luminosity function observed 
in the Milky Way and in other late type spirals. Even though our completeness 
limit is 5 x 10 erg s , in low density regions we are able to detect sources five 
times fainter than this, corresponding to the faintest possible Hn region. The 8 
and 24 luminosities within the central 5 kpc of M 33 are comparable and of or- 
der 4x 10 28 erg s^ 1 Hz" 1 (vL v (8) = 1. 5 xlO 42 and uL v (24) = 4.4x 10 41 erg s' 1 ). 
We estimate the total IR emission in the same region of M 33 to be 10 9 L . The 
discrete sources account for about one third of the 24 [im emission while the rest 
is diffuse. From the IR emission, we derive a star formation rate for the inner 
disk equal to 0.2 M Q yr _1 . 
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